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TWO PHASE CHOKE FLOW IN TUBES WITH VERY LARGE L/D 
by R. C. Hendricks and R. J. Simoneau 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

ABSTRACT 

Two phase and gaseous choked flow data for fluid nitrogen were obtained for a 
test section which was a long constant area duct of 16 200 L/D with a diverging 
diffuser attached to the exit. Flow rate data were taken along five isotherms (re- 
duced temperature of 0. 81, 0.96, 1.06, 1.12, and 2. 34) for reduced pressures 
to 3. The flow rate data were mapped in the usual manner using stagnation condi- 
tions at tlie inlet mixing chamber upstream of the entrance length. The results are 
predictable by a two-phase homogeneous equilibrium choking flow model which in- 
cludes wall friction. A simplified theory v^ich in essence decouples the long tube 
region from the high acceleration choking region also appears to predict the data 
reasonably well, but about 15 percent low. 
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TWO PHASE CHOKED FLOW IN TUBES WITH VERY LARGE L/D 
by R. C. Hendricks and R. J, Simoneau 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

INTRODUCTION 

Currently, the shuttle engine turbopump is required to boost propellant pres- 
sures to 30 MPa with proposed second generation engines requiring propellants to 
be delivered at pressures to 50 MPa. The problem of fluid leaking past the sealing 
surfaces in rotating machinery is compoimded with cryogenics, hi^ pressure, large 
temperature gradients, very high speeds of rotation, and static seal requirements. 

At lower pressure and rotation speeds, self energizing pumping seals with very 
close clearances have been successfully employed in a variety of sealing applications 
(Zuk et al. [^) ). These seals frequently have very large length to hydraulic diameter 
ratio (L/D) passages. At the proposed operating pressures, design innovations to 
minimize losses are required, but choked flow data and models to make such calcu- 
lations are lacking, A similar problem occurs in very long cryogenic transmission 
lines. 

Most two phase choked flow data reported in the literature are from experimen- 
tal devices with low L/D sections with little attention given to large L/D lines. 

The literature has been well surveyed in referenees 2 to 4. In previous experi- 
ments, the authors have studied two phase choked flows in a variety of geometries 
[ 1 including the orifice as a limiting case. References 12 to 16 have considered 

short tubes; in general for L/D < 3 short tubes behave much like the orifice and the 
effects of friction may be neglected. For larger L/D the effects of friction become 
increasingly important in determining the limiting mass flow rate and pressure drop. 
The question of two phase choked flows in large L/D tubes has, to the authors know- 
ledge, never been resolved. Toward this end, experimental two phase choked flow 
data for fluid nitrogen in a tube of 16 200 L/D are presented herein. 
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These results should have several applications including; aerospace, aeronau- 
tical and stationary engines where hi^er pressure components will be used to 
achieve hi^er efficiencies; high operating pressure evaporators, liquifiers and con- 
densers associated with the cryogenic and petrochemical industries and geothermal 
power production pipe lines for transmission of energy/power. In addition the re- 
sults should be an aid in defining the ultimate cooling capacity of a fluid in a heat ex- 
change device; and an aid in defining the nature of metastability and reasonable 
boundaries for metastable operation, 

DESCRIPTION OF THE APPARATUS AND PROCEDURE 

The flow system (fig. 1) is essentially that of reference 17 but modified for the 
present test. The test section was a coil which normally served as a heat shield for 
other test sections. By a fairly easy rearrangement of the plumbing tlie flow could 
be diverted through the coil. 

The 16 200 L/D test apparatus was made by winding a coil of 54 turns of 
0. 64 cm (0. 25 in. ) o. d. by 0. 48 cm (0. 19 in. ) i. d. copper tubing on a 45. 7 cm 
(18 in. ) diameter drum, yielding a 78. 3 meter tube length. Coil spacing is main- 
tained at 0. 64 cm by three bakelite strips which also support the coil. The exterior 
of the coil was covered with H cm of multilayer insulation and as shown in figure 1, 
the entire apparatus was located in an evacuated environment. 

The pressure taps were fabricated from 0, 32 cm (l/8 in. ) diameter tubing 
silver soldered to the tube along the inner surface of the coil at 10 coil intervals ex- 
cept for the last tap. The axial distance between static pressure measurements are 
given in Table I in terms of L/D with the first location on the tube taken as zero. 

In the manner suggested by the work of Henry [ | two diffusers were fabricated 

and attached to the end of the long tube, one with a 7° half angle divergence section 
and the other a 3^® half angle divergence section, see figure 2. It has been ob- 
served [ | that exit pressure measurement can be made more accurately if the exit 

is not abrupt. To access the effect of static tap diameter on the flow, two static taps 
were placed in the 7° half angle diffuser; a 0.04 cm (l/64 in. ) diameter hole was 
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located 0. 64 cm (| in. ) from the diverger.ce plane and a 0. 08 cm (1/16 in. ) diameter 
hole located at 0. 32 cm (l/8 in. ) from the divergence plane. No adverse effects due 
to the pressure taps were noted. The exit pressure tap, e in Table I, was less than 
1 diameter from the divergence plane. 

A static pressure was measured approximately 8 tube diameters downstream 
from the choking plane in the 1. 6 cm diameter transition tube at the exit of the dif- 
fuser leading to the back pressure control valve. This was designated back pressure. 

Fluici tenq>eratures were monitored at four positions with the following instru- 
ments; platinum thermometer in the mixing chamber, in the line open ball thermo- 
couples at approximately 75 L/D from the mixing chamber (inlet), and approxi- 
mately 50 L/D from the nozzle (outlet) and an (^en ball thermocouple located in a 
well in the diffuser body. 

The tank w'as filled with liquid nitrogen and pressurized. The normal exhaust 
valve was closed in order to divert the flow into the coil. The flow control valve 
was wide open, and the back pressure valve was opened; the flow rate throu^ the 
coil was measured at the venturi and the orifice. Adjustments of the back pressure 
valve caused changes in the static pressure within the transition tube without changes 
in the exit plane pressure or pressure profile which assured that the coil was choked. 

The inlet temperature to the coil could be increased using the joule heated tube 
(central test apparatus) which served as a preheater for these tests. 

ANALYSIS 

Consider a tube of length L and diameter D; 
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The tube has potentially three flow regions; 

Entrance (0-1): Since we are dealing with a single phase incompressible, or 
nearly incompressible, fluid at the entrance, we can simply assign a standard head 
loss for the entrance configuration. Then using the friction factor we can extend the 
tube length appropriately. 

Single phase (1-2): This region is normally dominated by friction but as the 
fluid gets closer to the thermodjmamic critical point the flow will become more com- 
pressible. In some cases points 2 and e will coincide. 

Two-phase (2-e): This region is normally dominated by the momentum of the 
expanding two-phase fluid, however, if vaporization occurs early enough in a long 
tube, wall friction may also be important in this region. 

The basic equation for one- dimensional flow with friction can be written 

dP + G^^+^dz = 0 (1) 

fv fv 2D 


where 


v(P) 


P ^ P- 




t m ' 


kx 


^ 1 -. 


sat 

xVg + (1 - x)kvjj 

P < P=«f 
sat 


for a given path (e. g. , isenthalpic) 
f = f(Re) = ff^\ = f(P) 


( 2 ) 


( 3 ) 


The friction factor f is the standard Fanning friction factor and is computed 
based on Reynolds number for a smooth tube. Since G is constant the friction 
factor is dependent only on the fluid viscosity and thus becomes a thermodynamic 
variable. In the two-phase region the question becomes one of what fluid is in 
contact with the wall. At least initially it is ejq>ected that there will be a vapor 
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annulus and liquid core ( ]. Later the flow becomes thoroughly dispersed and the 
overwhelming majority of the wall area will be in contact with vapor. Thus in this 
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analysis the two- phase wall friction factor will be based on the vapor viscosity. 

The flow is assumed to be homogeneous (k = 1) and in thermodynamic equilibrium. 
Reference 3 suggests that this assumption becomes more valid as tube length in- 
creases, which is the case herein. The thermodynamic path is assumed to be isen- 
thalpic. Since G is constant, equation (1) can be integrated and rearranged in the 
form 

P 


g2 = 


•Pn 


— dv 
fv 


j. ^ dp 


( 4 ) 


fv dP 


2D 


0 


To determine choking, the derivative dG/dP is examined at the exit plane. Differ- 
entiating equation (4) 3 delds 


2G~ 
dP ' 


J-^dP 

fv dP 


4- -5- = - ■ — — + -i- - ! — 

2D \ fv dP 2D dpj > fv 


( 5 ) 
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The choking condition requires dG/dP)^ = 0. Thus 

0 


max 


dP 2D dP 


(6) 


It is known that at the exit plane dP/dz becomes very large, approaching infinity 
in the ideal case. Thus at the exit dz/dP can be neglected relative to dv/dP with 
the familiar result 

-1 


g: 


max 


dv 

dP 


0 ) 


By substituting L' for z in equation (4) and iterating between equations (4) and (7) 
the exit pressure, P„, and the choked maximum flow rate, G^„„, can be obtained. 

0 TriuJC 

While in the present report the point of vaporization is taken as the equilibrium 
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saturation pressure, a nonequilibrium vaporization point could be defined. 

Once the maximum mass flux is determined, equation (4) can be used to solve 
for the pressure distribution along the tube. 

Special Case 

In many cases, especially in very long tubes and vv^en the initial conditions are 
well away from saturation, some simplifying assumptions can be applied, which 
allow simple estimates of the maximum flow rates. 

1. The two- phase region, dominated by momentum, is very short and friction 
can be neglected. 

2. The single-phase region, dominated by friction, is very long and relatively 
incompressible. Thus specific volume and friction factor can be taken as constant 
and the length Lq _2 can be taken as L'. 

Thus equation (4) can be simplified to 


pi- 

max 


2(Pq - Psat) 

vf(L'/D) 


( 8 ) 


Clearly equation (8) is the standard single phase friction flow equation found in any 
fluid mechanics text and it is not necessary to use the above analysis to arrive at it. 
It is recommended herein that the proper approach is the simultaneous solution of 
equations (4) and (7) with real fluid properties, such as available in reference 19. 

On the other hand it is reassuring that equations (4) and (7) reduce t*" a well-known 
form and that equation (8) can be used to compute flow rate in a fair number of 
cases. The exercise in obtaining the simplification also shows regimes where it 
might be applicable. 


RESULTS 


Data for fluid nitrogen were acquired and are presented in Table II. The five 
isotherms are presented as figure 3. In figure 3, the mass flux G and the reduced 
mass flux (0^. = G/G*) are given as a function of reduced pressure Pq/P^ 

where Pq is the pressure measured at the mixing chamber. The five isotherms 
(0. 81, 0. 96, 1. 06, 1. 12, and 2. 34) are also based on the mi.xing chamber tempera- 
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ture = Tq/T^. The data appear to be repre 3nted quite well using these mixing 
chamber parameters even thou^ the choked interface is over 16 200 L/D from the 
mixing chamber. 

The mass flux-pressure map of figure 3 resembles the shape and trends of the 
same map for a nozzle without an entry length at the same reduced pressure; how- 
ever, the value of _ is approximately 0. 1 that for the nozzle without the 

16 200 L/D entry section. 

The analysis tends to underpredict the flow at the high pressure end by about 
12 percent, while at the lower pressures the data and analysis are quite close. The 
program was designed for computing isenthalpic below the critical enthalpy; how- 
ever, some of the data at higher temperatures and lower pressures are above criti- 
cal enthalpy, thus only a partial check can be made. The difference in the trend is 
somewhat disturbing. It could be in the data. The data plotted were metered with 
the orifice downstream of the test section. The upstream venturi was about 8 per- 
cent lower at hi^ pressure but at low pressure was unreliable. Since the orifice 
showed a consistent trend throughout the authors stuck with it. It should be pointed 
out, however, that in these long tubes the flows were about one-tenth the normal flow 
for the test rig and thus the usual high level of flow accuracy was not present. It also 
should be pointed out that even thou^ there exists a 10 to 15 percent discrepancy in 
measured versus calculated flow the one- dimensional homogeneous equilibrium 
analysis with friction is closely predicting flows that are an order of magnitude below 
nozzle flows. 

Such a reduction is borne out by the pressure - L/D profiles of figure 4, where 
the pressure at 100 L/D ahead of tlic nozzle is approximately 0, 1 Pq. For the sake 
of clarity only six profiles were plotted; however, these are typical of the data pre- 
sented in figure 3. They are tabulated in Table III, Curiously the analysis does an 
excellent job of predicting the pressure profiles over the whole range of the experi- 
ment. 
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Another interesting trend appears in the slope of the pressure profiles depending 
on whether or not Tq is above or below T^, (T^ = 126. 3 K), Such is associated 
with the friction- momentum pressure losses in the coil as « 1. Actually a closer 
examination shows the profiles tend to become more curved as Hq -► All of the 
data are for Hq < however, the two dashed curves are for Hq very close to H^. 
The calculated curve for Tq = 141. 9 K actually crosses all the other curves and so 
do the data. At low temp>erature (and low enthalpy) the calculated profiles are linear 
in the single phase region, however, the data show a sli^t concave curvature. 

Figure 5 is a plot of the simple equation (8). The more exact analysis is super- 
imposed. As can be seen the equation of simple single phase friction does a good 
job, although it underpredicts fuilher. It can be used with confidence for quick cal- 
culations. Not too surprisingly, the simple friction approach is best at low temper- 
ature where the fluid is hi^ly incompressible and the saturation pressure is low. 

SUMMARY 

Two phase choked flow data for fluid nitrogen flowing throu^ a 16 2000 L/D 
coll were taken along five isotherms (0. 81, 0.96, 1.06, 1. 12, and 2. 34). These 
data can be represented by the same parameters associated with two phase choked 
flows without a large entry length. 

Theoretical calculations indicate that a homogeneous equilibrium model, repre- 
senting an integral of friction and momentum predicts the flow data with reasonable 
accuracy and the pressure profiles quite well. A simple equation based on the fric- 
tion pressure drop in the single phase region also predicts the flow quite well. 

NOI'ATION 

D diameter, c:m 

f friction factor (Darcy) 

9 

G flow rate, g/rm"'-s 

G* flow normalizing parameter (6010 g/cm^-s for nitrogen),^ ^c^c'^^c 
H enthalpy, J/g-K 
k two-phase slip ratio, uVu^ 
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L length, cm 

L' length, including adjustment for entrance losses, cm 

P pressure, MPa 

Q 

R gas constant, MPacm /g-K 

Re Reynolds number, GD/ju 

T temperature, K 

U velocity, cm/ s 

V specific volume, cm /g 

X quality 

Z compressibility, P/pRT 

z axial distance along tube, cm 

H viscosity, g/cm-s 

Q 

p density, g/cm 

Subscripts 

c thermodynamic critical point 

e tube exit 

g saturated gas (vapor) 

H saturated liquid 

max maximum value 

r reduced by corresponding states parameter 

sat saturation value 

0 stagnation (tube entrance) 

1 point on the tube 

2 point of two phase in the tube 

Superscripts 

* reference condition 

“ average value 
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Table n. Data Summary 
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Table ni. Selected Pressure Profiles 
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Figure 4 - Experi'^ental and calculated a>lal pressure profiles lor the 16 200 UD test coil. 











